Transcriptomic Data Analysis Using the Galaxy Platform: Coffee (Coffea arabica L.) Flowers as Example.

Coffee, an important agricultural product for tropical producing countries, is facing challenges due to climate change, including periods of drought, irregular rain distribution, and high temperatures. These changes result in plant water stress, leading to significant losses in coffee productivity and quality. Understanding the processes that affect coffee flowering is crucial for improving productivity and quality. In this chapter, we describe a protocol for transcriptome analysis using available Internet software, mainly in the Galaxy Platform, using RNA-Seq data from flowers collected from different parts of the coffee tree. The methods presented in this chapter provide a comprehensive protocol for transcriptome analysis of differentially expressed genes from flowers of coffee plant. This knowledge can be utilized in coffee genetic improvement programs, particularly in the selection of cultivars that are tolerant to water deficit.

Single-cell insights: pioneering an integrated atlas of chromatin accessibility and transcriptomic landscapes in diabetic cardiomyopathy.

BACKGROUND: Diabetic cardiomyopathy (DCM) poses a growing health threat, elevating heart failure risk in diabetic individuals. Understanding DCM is crucial, with fibroblasts and endothelial cells playing pivotal roles in driving myocardial fibrosis and contributing to cardiac dysfunction. Advances in Multimodal single-cell profiling, such as scRNA-seq and scATAC-seq, provide deeper insights into DCM's unique cell states and molecular landscape for targeted therapeutic interventions. METHODS: Single-cell RNA and ATAC data from 10x Multiome libraries were processed using Cell Ranger ARC v2.0.1. Gene expression and ATAC data underwent Seurat and Signac filtration. Differential gene expression and accessible chromatin regions were identified. Transcription factor activity was estimated with chromVAR, and Cis-coaccessibility networks were calculated using Cicero. Coaccessibility connections were compared to the GeneHancer database. Gene Ontology analysis, biological process scoring, cell-cell communication analysis, and gene-motif correlation was performed to reveal intricate molecular changes. Immunofluorescent staining utilized various antibodies on paraffin-embedded tissues to verify the findings. RESULTS: This study integrated scRNA-seq and scATAC-seq data obtained from hearts of WT and DCM mice, elucidating molecular changes at the single-cell level throughout the diabetic cardiomyopathy progression. Robust and accurate clustering analysis of the integrated data revealed altered cell proportions, showcasing decreased endothelial cells and macrophages, coupled with increased fibroblasts and myocardial cells in the DCM group, indicating enhanced fibrosis and endothelial damage. Chromatin accessibility analysis unveiled unique patterns in cell types, with heightened transcriptional activity in myocardial cells. Subpopulation analysis highlighted distinct changes in cardiomyocytes and fibroblasts, emphasizing pathways related to fatty acid metabolism and cardiac contraction. Fibroblast-centered communication analysis identified interactions with endothelial cells, implicating VEGF receptors. Endothelial cell subpopulations exhibited altered gene expressions, emphasizing contraction and growth-related pathways. Candidate regulators, including Tcf21, Arnt, Stat5a, and Stat5b, were identified, suggesting their pivotal roles in DCM development. Immunofluorescence staining validated marker genes of cell subpopulations, confirming PDK4, PPARγ and Tpm1 as markers for metabolic pattern-altered cardiomyocytes, activated fibroblasts and endothelial cells with compromised proliferation. CONCLUSION: Our integrated scRNA-seq and scATAC-seq analysis unveils intricate cell states and molecular alterations in diabetic cardiomyopathy. Identified cell type-specific changes, transcription factors, and marker genes offer valuable insights. The study sheds light on potential therapeutic targets for DCM.

Integration of single-cell RNA-seq and bulk RNA-seq to construct liver hepatocellular carcinoma stem cell signatures to explore their impact on patient prognosis and treatment.

BACKGROUND: Liver hepatocellular carcinoma (LIHC) is a prevalent form of primary liver cancer. Research has demonstrated the contribution of tumor stem cells in facilitating tumor recurrence, metastasis, and treatment resistance. Despite this, there remains a lack of established cancer stem cells (CSCs)-associated genes signatures for effectively predicting the prognosis and guiding the treatment strategies for patients diagnosed with LIHC. METHODS: The single-cell RNA sequencing (scRNA-seq) and bulk RNA transcriptome data were obtained based on public datasets and computerized firstly using CytoTRACE package and One Class Linear Regression (OCLR) algorithm to evaluate stemness level, respectively. Then, we explored the association of stemness indicators (CytoTRACE score and stemness index, mRNAsi) with survival outcomes and clinical characteristics by combining clinical information and survival analyses. Subsequently, weighted co-expression network analysis (WGCNA) and Cox were applied to assess mRNAsi-related genes in bulk LIHC data and construct a prognostic model for LIHC patients. Single-sample gene-set enrichment analysis (ssGSEA), Cell-type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORT) and Tumor Immune Estimation Resource (TIMER) analysis were employed for immune infiltration assessment. Finally, the potential immunotherapeutic response was predicted by the Tumor Immune Dysfunction and Exclusion (TIDE), and the tumor mutation burden (TMB). Additionally, pRRophetic package was applied to evaluate the sensitivity of high and low-risk groups to common chemotherapeutic drugs. RESULTS: A total of four genes (including STIP1, H2AFZ, BRIX1, and TUBB) associated with stemness score (CytoTRACE score and mRNAsi) were identified and constructed a risk model that could predict prognosis in LIHC patients. It was observed that high stemness cells occurred predominantly in the late stages of LIHC and that poor overall survival in LIHC patients was also associated with high mRNAsi scores. In addition, pathway analysis confirmed the biological uniqueness of the two risk groups. Personalized treatment predictions suggest that patients with a low risk benefited more from immunotherapy, while those with a high risk group may be conducive to chemotherapeutic drugs. CONCLUSION: The current study developed a novel prognostic risk signature with genes related to CSCs, which provides novel ideas for the diagnosis, prognosis and treatment of LIHC.

Systematic mapping of TF-mediated cell fate changes by a pooled induction coupled with scRNA-seq and multi-omics approaches.

Transcriptional regulation controls cellular functions through interactions between transcription factors (TFs) and their chromosomal targets. However, understanding the fate conversion potential of multiple TFs in an inducible manner remains limited. Here, we introduce iTF-seq as a method for identifying individual TFs that can alter cell fate toward specific lineages at a single-cell level. iTF-seq enables time course monitoring of transcriptome changes, and with biotinylated individual TFs, it provides a multi-omics approach to understanding the mechanisms behind TF-mediated cell fate changes. Our iTF-seq study in mouse embryonic stem cells identified multiple TFs that trigger rapid transcriptome changes indicative of differentiation within a day of induction. Moreover, cells expressing these potent TFs often show a slower cell cycle and increased cell death. Further analysis using bioChIP-seq revealed that GCM1 and OTX2 act as pioneer factors and activators by increasing gene accessibility and activating the expression of lineage specification genes during cell fate conversion. iTF-seq has utility in both mapping cell fate conversion and understanding cell fate conversion mechanisms.

RNA-seq analysis-based study on the effects of gestational diabetes mellitus on macrosomia.

Background: Both the mother and the infant are negatively impacted by macrosomia. Macrosomia is three times as common in hyperglycemic mothers as in normal mothers. This study sought to determine why hyperglycemic mothers experienced higher macrosomia. Methods: Hematoxylin and Eosin staining was used to detect the placental structure of normal mother(NN), mothers who gave birth to macrosomia(NM), and mothers who gave birth to macrosomia and had hyperglycemia (DM). The gene expressions of different groups were detected by RNA-seq. The differentially expressed genes (DEGs) were screened with DESeq2 R software and verified by qRT-PCR. The STRING database was used to build protein-protein interaction networks of DEGs. The Cytoscape was used to screen the Hub genes of the different group. Results: The NN group's placental weight differed significantly from that of the other groups. The structure of NN group's placenta is different from that of the other group, too. 614 and 3207 DEGs of NM and DM, respectively, were examined in comparison to the NN group. Additionally, 394 DEGs of DM were examined in comparison to NM. qRT-PCR verified the results of RNA-seq. Nucleolar stress appears to be an important factor in macrosomia, according on the results of KEGG and GO analyses. The results revealed 74 overlapped DEGs that acted as links between hyperglycemia and macrosomia, and 10 of these, known as Hub genes, were key players in this process. Additionally, this analysis believes that due of their close connections, non-overlapping Hubs shouldn't be discounted. Conclusion: In diabetic mother, ten Hub genes (RPL36, RPS29, RPL8 and so on) are key factors in the increased macrosomia in hyperglycemia. Hyperglycemia and macrosomia are linked by 74 overlapping DEGs. Additionally, this approach contends that non-overlapping Hubs shouldn't be ignored because of their tight relationships.

Genome-wide annotation of transcript boundaries using bacterial Rend-seq datasets.

Accurate annotation to single-nucleotide resolution of the transcribed regions in genomes is key to optimally analyse RNA-seq data, understand regulatory events and for the design of experiments. However, currently most genome annotations provided by GenBank generally lack information about untranslated regions. Additionally, information regarding genomic locations of non-coding RNAs, such as sRNAs, or anti-sense RNAs is frequently missing. To provide such information, diverse RNA-seq technologies, such as Rend-seq, have been developed and applied to many bacterial species. However, incorporating this vast amount of information into annotation files has been limited and is bioinformatically challenging, resulting in UTRs and other non-coding elements being overlooked or misrepresented. To overcome this problem, we present pyRAP (python Rend-seq Annotation Pipeline), a software package that analyses Rend-seq datasets to accurately resolve transcript boundaries genome-wide. We report the use of pyRAP to find novel transcripts, transcript isoforms, and RNase-dependent sRNA processing events. In Bacillus subtilis we uncovered 63 novel transcripts and provide genomic coordinates with single-nucleotide resolution for 2218 5'UTRs, 1864 3'UTRs and 161 non-coding RNAs. In Escherichia coli, we report 117 novel transcripts, 2429 5'UTRs, 1619 3'UTRs and 91 non-coding RNAs, and in Staphylococcus aureus, 16 novel transcripts, 664 5'UTRs, 696 3'UTRs, and 81 non-coding RNAs. Finally, we use pyRAP to produce updated annotation files for B. subtilis 168, E. coli K-12 MG1655, and S. aureus 8325 for use in the wider microbial genomics research community.

RNA Isolation from Ctenophores.

RNA-seq or transcriptome analysis of individual cells and small cell populations is essential for virtually any biomedical field. Here, we examine and discuss the different methods of RNA isolation specific to ctenophores. We present a convenient, inexpensive, and reproducible protocol for RNA-seq libraries that are designed for low quantities of samples. We demonstrated these methods on early (one, two, four, eight cells) embryonic and developmental stages, tissues, and even a single aboral organ from the ctenophore Pleurobrachia bachei and other ctenophore species (e.g., Mnemiopsis, Bolinopsis, and Beroe).

Ocean to Tree: Leveraging Single-Molecule RNA-Seq to Repair Genome Gene Models and Improve Phylogenomic Analysis of Gene and Species Evolution.

Understanding gene evolution across genomes and organisms, including ctenophores, can provide unexpected biological insights. It enables powerful integrative approaches that leverage sequence diversity to advance biomedicine. Sequencing and bioinformatic tools can be inexpensive and user-friendly, but numerous options and coding can intimidate new users. Distinct challenges exist in working with data from diverse species but may go unrecognized by researchers accustomed to gold-standard genomes. Here, we provide a high-level workflow and detailed pipeline to enable animal collection, single-molecule sequencing, and phylogenomic analysis of gene and species evolution. As a demonstration, we focus on (1) PacBio RNA-seq of the genome-sequenced ctenophore Mnemiopsis leidyi, (2) diversity and evolution of the mechanosensitive ion channel Piezo in genetic models and basal-branching animals, and (3) associated challenges and solutions to working with diverse species and genomes, including gene model updating and repair using single-molecule RNA-seq. We provide a Python Jupyter Notebook version of our pipeline (GitHub Repository: Ctenophore-Ocean-To-Tree-2023 https://github.com/000generic/Ctenophore-Ocean-To-Tree-2023 ) that can be run for free in the Google Colab cloud to replicate our findings or modified for specific or greater use. Our protocol enables users to design new sequencing projects in ctenophores, marine invertebrates, or other novel organisms. It provides a simple, comprehensive platform that can ease new user entry into running their evolutionary sequence analyses.

Immune microniches shape intestinal Treg function.

The intestinal immune system is highly adapted to maintaining tolerance to the commensal microbiota and self-antigens while defending against invading pathogens1,2. Recognizing how the diverse network of local cells establish homeostasis and maintains it in the complex immune environment of the gut is critical to understanding how tolerance can be re-established following dysfunction, such as in inflammatory disorders. Although cell and molecular interactions that control T regulatory (Treg) cell development and function have been identified3,4, less is known about the cellular neighbourhoods and spatial compartmentalization that shapes microorganism-reactive Treg cell function. Here we used in vivo live imaging, photo-activation-guided single-cell RNA sequencing5-7 and spatial transcriptomics to follow the natural history of T cells that are reactive towards Helicobacter hepaticus through space and time in the settings of tolerance and inflammation. Although antigen stimulation can occur anywhere in the tissue, the lamina propria-but not embedded lymphoid aggregates-is the key microniche that supports effector Treg (eTreg) cell function. eTreg cells are stable once their niche is established; however, unleashing inflammation breaks down compartmentalization, leading to dominance of CD103+SIRPα+ dendritic cells in the lamina propria. We identify and validate the putative tolerogenic interaction between CD206+ macrophages and eTreg cells in the lamina propria and identify receptor-ligand pairs that are likely to govern the interaction. Our results reveal a spatial mechanism of tolerance in the lamina propria and demonstrate how knowledge of local interactions may contribute to the next generation of tolerance-inducing therapies.

RNA-Seq analysis reveals the mechanism in response to cold stress of peach cv 'Dingjiaba'.

BACKGROUND: 'Dingjiaba' is an important Prunus persica cultivar (cv) mainly grown in the Hexi corridor in northwest China, which has an inherited strong cold tolerance. OBJECTIVE: To compare the transcriptome and physiology data of leaves of cvs 'Dingjiaba' (D) and 'Kanoiwa' (K) following cold treatment at different time periods, in order to gain new insights into the mechanisms of cold adaptation in 'Dingjiaba'. MATERIALS AND METHODS: We analyzed the transcriptomic and physiological data of leaves of D and K cvs exposed to 0 h (D0/K0), 2 h (D2/K2), 6 h (D6/K6) and 12 h (D12/K12) cold stress. RESULTS: Low temperature stress caused membrane damage and led to increased rate of electrolyte leakage and increased MDA content. Cold stress induced the accumulation of soluble sugars, soluble proteins and proline in leaves of both cvs, with a lower increase in K compared to D. Transcriptome analysis identified 4,631, 5,069, 5,662 and 3,886 differentially expressed genes (DEGs) between D0 and K0, D2 and K2, D6 and K6 and D12 and K12, respectively. The differentially expressed genes significantly enriched in metabolic pathways and biosynthesis of secondary metabolites. We further validated the reliability of sequencing data of the RNA-Seq with Real-Time Quantitative PCR, which suggested that the expression trend of the RNA-Seq were same as RT-PCR. CONCLUSIONS: These results provide novel insights into a series of molecular mechanisms underlying physiological metabolism and defense. https://doi.org/10.54680/fr24210110312.

A concerted increase in readthrough and intron retention drives transposon expression during aging and senescence.

Aging and senescence are characterized by pervasive transcriptional dysfunction, including increased expression of transposons and introns. Our aim was to elucidate mechanisms behind this increased expression. Most transposons are found within genes and introns, with a large minority being close to genes. This raises the possibility that transcriptional readthrough and intron retention are responsible for age-related changes in transposon expression rather than expression of autonomous transposons. To test this, we compiled public RNA-seq datasets from aged human fibroblasts, replicative and drug-induced senescence in human cells, and RNA-seq from aging mice and senescent mouse cells. Indeed, our reanalysis revealed a correlation between transposons expression, intron retention, and transcriptional readthrough across samples and within samples. Both intron retention and readthrough increased with aging or cellular senescence and these transcriptional defects were more pronounced in human samples as compared to those of mice. In support of a causal connection between readthrough and transposon expression, analysis of models showing induced transcriptional readthrough confirmed that they also show elevated transposon expression. Taken together, our data suggest that elevated transposon reads during aging seen in various RNA-seq dataset are concomitant with multiple transcriptional defects. Intron retention and transcriptional readthrough are the most likely explanation for the expression of transposable elements that lack a functional promoter.

Shaoxia: a web-based interactive analysis platform for single cell RNA sequencing data.

BACKGROUND: In recent years, Single-cell RNA sequencing (scRNA-seq) is increasingly accessible to researchers of many fields. However, interpreting its data demands proficiency in multiple programming languages and bioinformatic skills, which limited researchers, without such expertise, exploring information from scRNA-seq data. Therefore, there is a tremendous need to develop easy-to-use software, covering all the aspects of scRNA-seq data analysis. RESULTS: We proposed a clear analysis framework for scRNA-seq data, which emphasized the fundamental and crucial roles of cell identity annotation, abstracting the analysis process into three stages: upstream analysis, cell annotation and downstream analysis. The framework can equip researchers with a comprehensive understanding of the analysis procedure and facilitate effective data interpretation. Leveraging the developed framework, we engineered Shaoxia, an analysis platform designed to democratize scRNA-seq analysis by accelerating processing through high-performance computing capabilities and offering a user-friendly interface accessible even to wet-lab researchers without programming expertise. CONCLUSION: Shaoxia stands as a powerful and user-friendly open-source software for automated scRNA-seq analysis, offering comprehensive functionality for streamlined functional genomics studies. Shaoxia is freely accessible at http://www.shaoxia.cloud , and its source code is publicly available at https://github.com/WiedenWei/shaoxia .

Preparation of Healthy Single Neuron or Astrocyte Suspension from Adult Mouse Brain for RNA-seq.

Single-cell RNA sequencing (scRNA-seq) has been widely applied in neuroscience research, enabling the investigation of cellular heterogeneity at the transcriptional level, the characterization of rare cell types, and the detailed analysis of the stochastic nature of gene expression. Isolation of single nerve cells in good health, especially from the adult rodent brain, is the most difficult and critical process for scRNA-seq. Here, we describe methods to optimize protease digestion of brain slices, which enable yield of millions of cells in good health from the adult brain.

The landscape of abiotic and biotic stress-responsive splice variants with deep RNA-seq datasets in hot pepper.

Alternative splicing (AS) is a widely observed phenomenon in eukaryotes that plays a critical role in development and stress responses. In plants, the large number of RNA-seq datasets in response to different environmental stressors can provide clues for identification of condition-specific and/or common AS variants for preferred agronomic traits. We report RNA-seq datasets (350.7 Gb) from Capsicum annuum inoculated with one of three bacteria, one virus, or one oomycete and obtained additional existing transcriptome datasets. In this study, we investigated the landscape of AS in response to environmental stressors, signaling molecules, and tissues from 425 total samples comprising 841.49 Gb. In addition, we identified genes that undergo AS under specific and shared stress conditions to obtain potential genes that may be involved in enhancing tolerance to stressors. We uncovered 1,642,007 AS events and identified 4,354 differential alternative splicing genes related to environmental stressors, tissues, and signaling molecules. This information and approach provide useful data for basic-research focused on enhancing tolerance to environmental stressors in hot pepper or establishing breeding programs.

scINRB: single-cell gene expression imputation with network regularization and bulk RNA-seq data.

Single-cell RNA sequencing (scRNA-seq) facilitates the study of cell type heterogeneity and the construction of cell atlas. However, due to its limitations, many genes may be detected to have zero expressions, i.e. dropout events, leading to bias in downstream analyses and hindering the identification and characterization of cell types and cell functions. Although many imputation methods have been developed, their performances are generally lower than expected across different kinds and dimensions of data and application scenarios. Therefore, developing an accurate and robust single-cell gene expression data imputation method is still essential. Considering to maintain the original cell-cell and gene-gene correlations and leverage bulk RNA sequencing (bulk RNA-seq) data information, we propose scINRB, a single-cell gene expression imputation method with network regularization and bulk RNA-seq data. scINRB adopts network-regularized non-negative matrix factorization to ensure that the imputed data maintains the cell-cell and gene-gene similarities and also approaches the gene average expression calculated from bulk RNA-seq data. To evaluate the performance, we test scINRB on simulated and experimental datasets and compare it with other commonly used imputation methods. The results show that scINRB recovers gene expression accurately even in the case of high dropout rates and dimensions, preserves cell-cell and gene-gene similarities and improves various downstream analyses including visualization, clustering and trajectory inference.

The integration of single-cell and bulk RNA-seq atlas reveals ERS-mediated acinar cell damage in acute pancreatitis.

BACKGROUND: Acute pancreatitis (AP) is a clinically common acute abdominal disease, whose pathogenesis remains unclear. The severe patients usually have multiple complications and lack specific drugs, leading to a high mortality and poor outcome. Acinar cells are recognized as the initial site of AP. However, there are no precise single-cell transcriptomic profiles to decipher the landscape of acinar cells during AP, which are the missing pieces of jigsaw we aimed to complete in this study. METHODS: A single-cell sequencing dataset was used to identify the cell types in pancreas of AP mice and to depict the transcriptomic maps in acinar cells. The pathways' activities were evaluated by gene sets enrichment analysis (GSEA) and single-cell gene sets variation analysis (GSVA). Pseudotime analysis was performed to describe the development trajectories of acinar cells. We also constructed the protein-protein interaction (PPI) network and identified the hub genes. Another independent single-cell sequencing dataset of pancreas samples from AP mice and a bulk RNA sequencing dataset of peripheral blood samples from AP patients were also analyzed. RESULTS: In this study, we identified genetic markers of each cell type in the pancreas of AP mice based on single-cell sequencing datasets and analyzed the transcription changes in acinar cells. We found that acinar cells featured acinar-ductal metaplasia (ADM), as well as increased endocytosis and vesicle transport activity during AP. Notably, the endoplasmic reticulum stress (ERS) and ER-associated degradation (ERAD) pathways activated by accumulation of unfolded/misfolded proteins in acinar cells could be pivotal for the development of AP. CONCLUSION: We deciphered the distinct roadmap of acinar cells in the early stage of AP at single-cell level. ERS and ERAD pathways are crucially important for acinar homeostasis and the pathogenesis of AP.

BEERS2: RNA-Seq simulation through high fidelity in silico modeling.

Simulation of RNA-seq reads is critical in the assessment, comparison, benchmarking and development of bioinformatics tools. Yet the field of RNA-seq simulators has progressed little in the last decade. To address this need we have developed BEERS2, which combines a flexible and highly configurable design with detailed simulation of the entire library preparation and sequencing pipeline. BEERS2 takes input transcripts (typically fully length messenger RNA transcripts with polyA tails) from either customizable input or from CAMPAREE simulated RNA samples. It produces realistic reads of these transcripts as FASTQ, SAM or BAM formats with the SAM or BAM formats containing the true alignment to the reference genome. It also produces true transcript-level quantification values. BEERS2 combines a flexible and highly configurable design with detailed simulation of the entire library preparation and sequencing pipeline and is designed to include the effects of polyA selection and RiboZero for ribosomal depletion, hexamer priming sequence biases, GC-content biases in polymerase chain reaction (PCR) amplification, barcode read errors and errors during PCR amplification. These characteristics combine to make BEERS2 the most complete simulation of RNA-seq to date. Finally, we demonstrate the use of BEERS2 by measuring the effect of several settings on the popular Salmon pseudoalignment algorithm.

Deciphering antifungal and antibiofilm mechanisms of isobavachalcone against Cryptococcus neoformans through RNA-seq and functional analyses.

Cryptococcus neoformans has been designated as critical fungal pathogens by the World Health Organization, mainly due to limited treatment options and the prevalence of antifungal resistance. Consequently, the utilization of novel antifungal agents is crucial for the effective treatment of C. neoformans infections. This study exposed that the minimum inhibitory concentration (MIC) of isobavachalcone (IBC) against C. neoformans H99 was 8 µg/mL, and IBC dispersed 48-h mature biofilms by affecting cell viability at 16 µg/mL. The antifungal efficacy of IBC was further validated through microscopic observations using specific dyes and in vitro assays, which confirmed the disruption of cell wall/membrane integrity. RNA-Seq analysis was employed to decipher the effect of IBC on the C. neoformans H99 transcriptomic profiles. Real-time quantitative reverse transcription PCR (RT-qPCR) analysis was performed to validate the transcriptomic data and identify the differentially expressed genes. The results showed that IBC exhibited various mechanisms to impede the growth, biofilm formation, and virulence of C. neoformans H99 by modulating multiple dysregulated pathways related to cell wall/membrane, drug resistance, apoptosis, and mitochondrial homeostasis. The transcriptomic findings were corroborated by the antioxidant analyses, antifungal drug sensitivity, molecular docking, capsule, and melanin assays. In vivo antifungal activity analysis demonstrated that IBC extended the lifespan of C. neoformans-infected Caenorhabditis elegans. Overall, the current study unveiled that IBC targeted multiple pathways simultaneously to inhibit growth significantly, biofilm formation, and virulence, as well as to disperse mature biofilms of C. neoformans H99 and induce cell death.

Genome-wide characterization of post-transcriptional processes related to wood formation in Dalbergia odorifera.

BACKGROUND: Alternative polyadenylation (APA), alternative splicing (AS), and long non-coding RNAs (lncRNAs) play regulatory roles in post-transcriptional processes in plants. However, little is known about their involvement in xylem development in Dalbergia odorifera, a valuable rosewood species with medicinal and commercial significance. We addressed this by conducting Isoform Sequencing (Iso-Seq) using PacBio's SMRT technology and combined it with RNA-seq analysis (RNA sequencing on Illumina platform) after collecting xylem samples from the transition zone and the sapwood of D. odorifera. RESULTS: We identified 14,938 full-length transcripts, including 9,830 novel isoforms, which has updated the D. odorifera genome annotation. Our analysis has revealed that 4,164 genes undergo APA, whereas 3,084 genes encounter AS. We have also annotated 118 lncRNAs. Furthermore, RNA-seq analysis identified 170 differential alternative splicing (DAS) events, 344 genes with differential APA site usage (DE-APA), and 6 differentially expressed lncRNAs in the transition zone when compared to the sapwood. AS, APA, and lncRNAs are differentially regulated during xylem development. Differentially expressed APA genes were enriched for terpenoid and flavonoid metabolism, indicating their role in the heartwood formation. Additionally, DE-APA genes were associated with cell wall biosynthesis and terpenoid metabolism, implying an APA's role in wood formation. A DAS gene (involved in chalcone accumulation) with a significantly greater inclusion of the last exon in the transition zone than in the sapwood was identified. We also found that differentially expressed lncRNAs targeted the genes related to terpene synthesis. CONCLUSIONS: This study enhances our understanding of the molecular regulatory mechanisms underlying wood formation in D. odorifera, and provides valuable genetic resources and insights for its molecular-assisted breeding.

Integration of single-cell RNA-seq and bulk RNA-seq data to construct and validate a cancer-associated fibroblast-related prognostic signature for patients with ovarian cancer.

BACKGROUND: To establish a prognostic risk profile for ovarian cancer (OC) patients based on cancer-associated fibroblasts (CAFs) and gain a comprehensive understanding of their role in OC progression, prognosis, and therapeutic efficacy. METHODS: Data on OC single-cell RNA sequencing (scRNA-seq) and total RNA-seq were collected from the GEO and TCGA databases. Seurat R program was used to analyze scRNA-seq data and identify CAFs clusters corresponding to CAFs markers. Differential expression analysis was performed on the TCGA dataset to identify prognostic genes. A CAF-associated risk signature was designed using Lasso regression and combined with clinicopathological variables to develop a nomogram. Functional enrichment and the immune landscape were also analyzed. RESULTS: Five CAFs clusters were identified in OC using scRNA-seq data, and 2 were significantly associated with OC prognosis. Seven genes were selected to develop a CAF-based risk signature, primarily associated with 28 pathways. The signature was a key independent predictor of OC prognosis and relevant in predicting the results of immunotherapy interventions. A novel nomogram combining CAF-based risk and disease stage was developed to predict OC prognosis. CONCLUSION: The study highlights the importance of CAFs in OC progression and suggests potential for innovative treatment strategies. A CAF-based risk signature provides a highly accurate prediction of the prognosis of OC patients, and the developed nomogram shows promising results in predicting the OC prognosis.